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Abstract.  We  describe  an  approximation  scheme  which  can  be  used  to  estimate  unknown  parameters 
f  in  moving  boundary  problems.  The  model  equations  we  consider  are  fairly  general  nonlinear 
diffusion/reaction  equations  of  one  spatial  variable.  Here  we  give  conditions  on  the  parameter  sets  and 
model  equations  under  which  we  can  prove  that  the  estimates  obtained  using  the  approximations  will 
converge  to  best-fit  parameters  for  the  original  model  equations.  We  conclude  with  a  numerical 
example. 


Moving  boundary  problems  appear  in  a  large  variety  of  applications.  For  example,  in  [4],  a 
nonlinear  diffusion  equation  with  moving  boundary  is  derived  to  model  the  effects  of  biofilm  growth  on 
the  adsorption  properties  of  carbon  particles  (activated  carbon  filters  are  widely  used  to  remove 
pollutants  from  water).  In  [7],  parabolic  partial  differential  equations  involving  moving  boundaries  are 
employed  to  model  the  elongation  of  the  acrosomal  process  (a  key  step  in  the  fertilization  of  an  egg). 
Finally,  the  one  dimensional,  one  phase  Stefan  problem  is  well  known  as  a  model  for  the  melting  of  ice 
(see,  e.g.,  [5]). 

There  are  many  approaches  to  solving  these  problems,  both  analytical  and  numerical  (see,  e.g., 
[5],  [1]).  We  are  interested  in  numerical  methods  for  the  parameter  estimation,  or  inverse  problem;  we 
collect  data  by  observing  a  process  which  we  assume  can  be  modeled  by  a  moving  boundary  problem. 
We  then  try  to  estimate  unknown  parameters  appearing  in  the  model  equations  by  minimizing  a  least- 
squares  fit-to-data  criterion. 

Our  approach  to  solving  the  parameter  estimation  problem  is  that  of  [3]  (see  also  [2]  and  the 
references  of  both).  The  infinite  dimensional  model  equations  are  replaced  by  an  approximating  system 
of  ordinary  differential  equations  which  have  attractive  computational  features.  We  estimate  the 
unknown  parameters  within  these  approximating  equations.  We  prove  that  the  estimates  obtained 
using  the  approximate  model  equations  converge  in  an  appropriate  sense  to  best-fit  parameters  for  the 
original  equations. 

Specifically,  suppose  the  process  in  question  can  be  modeled  by  a  moving  boundary  problem 
having  the  following  form: 

ut  =  (^)(t,x)  ux  +  Y(t,x)u)  +  p(u)  +  f(t,x)  0  <  x  <  s(t),  0  <  t  <  T; 

u(0,x)  =  u0(x)  0  <  x  <  s0. 

(1)  «n(^  ux  +  rU)|  -  a12u(t,0)  =  g(t), 

v  /!x=0 

a21(<5  ux  +  Vu)|  +  ar22u(t,s(t))  =  h(l)  0  <  t  <  T; 

v  /lx=sCt) 

lj§  =  ?(s,u;7)  0  <  t  <  T  ;  s(0)  =  s0 

where  one  of  or  a12,  and  one  of  a21  or  a22  may  be  zero,  and  the  others  are 

positive;  the  functional  T  may  have  a  variety  of  forms,  to  be  discussed  in  more  detail  below,  and  7 
represents  a  parameter  which  is,  in  general,  a  function  of  (t,x).  We  assume  u0,  f,  g  and  h  are  known; 
the  unknowns  in  the  above  system  might  be  any  of  Y,  p,  Oy,  7,  or  s0,  denoted  in  the  following  by 
the  vector  q.  Let  us  assume  here  that  the  observations  are  in  the  form  of  point  evaluations,  i.e.,  we 
have  data  {  u j  j  }  corresponding  to  u(t|,Xj),  where  (u,s)  is  a  solution  of  (1).  Then  we  can  formulate 
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the  parameter  estimation  problem  as: 

m  n  2 

C?)  Minimize  J(q)  =  £  £  u,  ,  —  u(tj,xj;q)  for  qgQ  subject  to  (u,s) 

j=ii=r  1 

a  solution  of  (1). 

Given  the  general  model  discussed  above,  we  require  further  that  the  following  hypotheses  be  satisfied 
in  order  to  prove  convergence  of  the  parameter  estimates. 

(HI)  We  assume  that  (1)  has  a  unique  solution  (u,s)  €  cf  [0,T];  H1(0,s(t))^  x  H1(0,T] 
with  u(t,x)  >  0  and  s(t)  >  s  >  0.  ' 

(H2)  We  assume,  for  any  7,  7  satisfies  ^(s.Uj)  —  7(s,u2)|  <  A(s) jux  —  u2(  and  ^ 

|9(s1,u)  —  «F(s2,u)|  <  /i(|u|00)|s1  ~  s2|  with  A  and  p.  both  continuous.  * 

(H3)  We  shall  search  for  the  parameters  “5,  V,  and  7  within  compact  subsets  of 

C([0,T]  x  (0,s(t)]),  with  the  additional  requirement  that  there  exists  a  positive 
constant  d  such  that  all  satisfy  >  d;  we  search  for  p  within  a 

compact  subset  of  C(0,u)  (u  is  an  a  priori  upper  bound  on  u)  with  the 
additional  requirement  that  there  is  a  constant  L  such  that  all  p  satisfy 
|p(^i)  —  p(0 2)!  <  L|0i  —  02|  f°r  ^  R;  all  Qy  and  s0  belong  to  compact  subsets 
of  R+. 

The  second  hypothesis  is  satisfied,  for  example,  by  the  model  equations  for  the 
activated  carbon  problem  of  [4],  in  which  5  involves  an  integral  over  the  spatial  variable  of  a  nonlinear 
function  of  u  (see  the  numerical  example  below).  This  hypothesis  is  also  satisfied  by  one  of  the  models 
of  [7]  for  elongation  of  the  acrosomal  process,  in  which  5  involves  u(t,s(t)).  It  is  not  satisfied,  however 
by  the  Stefan  problem,  in  which  9  involves  ux(t,s(t)). 

In  order  to  derive  our  approximating  equations  and  to  facilitate  convergence  arguments,  we 
transform  the  original  problem  to  one  of  fixed  extent,  and  then,  as  in  [3],  we  rewrite  the  equation  in 
variational  form.  The  transformation  is  accomplished  by  setting  y  =  x/s(t),  letting  U(t,y)  =  u(t,x), 
U0(y)  =  u0(x),  F(t,y)  =  f(t,x),  D(t,y)  =  ^(t,x),  V(t,y)  =  f(t,x),  and  T(t,y)  =  7(t,x).  Let  us 
assume  for  definiteness  that  all  boundary  parameters  are  nonzero,  and  =  »21  =  1.  We  define 

X  =  H1(0,1)  and  let  ( • ,  • )  denote  the  inner  product  in  H°(0,1).  We  then  replace  equations  (1)  by: 

(ut,  0)  =  -4j  (DUy,  0y>  -  |  (V  U,  0y>  +  |  (y  Uy,  0}  +  (p(U),  0)  +  (F,  0) 

-  5  (g(t)  +  o,12U(t,0))  0(0)  +  l  (h(t)  -  o22U(t,l))  0(1) 

(2)  for  all  0  6  X,  t  €  [0,T] 

U(0,y)  =  U0(y) 

s  =  ?(s,  U;T)  0  <  t  <  Tj  s(0)  =  s0 

We  define  an  approximating  subspace  XN  C  X  as  the  span  of  the  set  of  cubic  B-splines  defined  on  a 
uniform  mesh  of  [0,1]  (see  [8]),  and  let  PN:  X  — >  XN  be  the  orthogonal  projection  in  the  H°  topology. 
Our  approximations  will  then  be  the  solution  (UN,  sN)  to  the  coupled  system: 

<u£\  0>  =  -pta  (DUy>  0y>  -  4j  (v  UN,  0y)  +  sJJ  (y  u£,  0)  +  <?(UN),  0)  +  (F,  0) 

^8  ^88 

-  4]  («(l)  +  al2UN(t.O))  0(0)  +  4q  (h(t)  -  a22UN(t,l))  0(1) 

8  8 

(2n)  for  all  0  e  XN,  t  €  [0,T] 

UN(0,y)  =  PNU0 

sN  =  «I(sN,UN;r)  0  <  t  <  T;  sN(0)  =  s0 
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Any  unknown  variable  coefficients  must  be  approximated  in  addition  to  the  state 
approximation  described  above.  Just  as  is  described  in  [2],  [3],  we  approximate  variable  parameters 
with  linear  splines,  replacing  the  parameter  set  Q  by  an  appropriate  finite  dimensional  approximating 
set  Qm.  The  corresponding  computationally  tractable  parameter  estimation  problem  is  given  by 


(cpN.M) 


Minimize  JN(q)  =  £  53 
j=li=l 


aiJ"UN(ti’^',q)! 


for  q  €  Qm  subject  to  (UN,sN) 


a  solution  of  (2^). 

The  system  of  equations  (2N)  represents  a  system  of  ordinary  differential  equations  which  can  be 
implemented  very  efficiently  on  the  computer  (the  details  of  the  numerical  aspects  will  appear  ^in^a 
forthcoming  manuscript;  a  limited  discussion  has  appeared  in  [6]).  The  optimization  problems  (9  '  ) 
are  solved  for  various  values  of  N  and  M,  resulting  in  a  sequence  of  best-fit  estimates.  As  in  [2],  [3]  the 
crucial  result  in  proving  that  solutions  of  (S’  '  )  converge  (subsequentially) 
to  a  solution  of  (S’)  is  that  ^UN(qN,M),  sN(qN,M)^  converges  (in  an  appropriate  sense)  to  ^U(q),  s(q)^ 
for  any  sequence  qN,M  which  converges  to  q  (in  an  appropriate  sense).  Such  a  convergence  statement 


can  be  proven  under  the  hypotheses  (HI)  -  (H3)  above,  with  the  addition  of: 


(H4)  There  exist  constants  s,  K  such  that  sN(t)  >  s  >  0  and  |sN(t)|  <  K  for  all  N 
and  t  €  [0,TJ. 


The  above  hypothesis  can  be  met  for  the  same  two  problems  mentioned  after  (H2).  Given  a  priori 
knowledge  about  solutions  of  (1),  one  can  modify  the  equation  for  sN  by  replacing  UN  with  a 
constrained  version,  UN;  this  constrained  approximate  state  is  constructed  in  such  a  way  that  both 
(114)  is  satisfied,  and  |u  -  UM|o0  <  ilJ  -  UN|OT  (thereby  ensuring  that  ?(sN,UN)  is  a  good 
approximation  to  <J(sN,UN)  ).  The  details  of  these  arguments  will  appear  elsewhere. 

We  present  the  results  of  a  test  problem,  motivated  by  the  activated  carbon  model  of  [4).  Our 
equations  are  (1)  with 

*  =  ]o  to  “  7(t))dx, 

V  =  0,  p  =  0,  au  =  1,  a12  =  0,  a21  =  1,  and  o22  =  1-  We  give  the  results  of  two  estimations, 
both  of  an  unknown  diffusion  coefficient.  In  both  cases,  our  data  is  taken  to  be  u(tj,0)  for  ten  values  of 
time.  To  generate  data,  we  have  first  chosen  functions  u  and  s,  and  a  "true”  diffusion  coefficient  1, 
and  then  determined  f,  g,  h,  7,  u0,  and  s0  so  that  (1)  holds.  We  assume  all  parameters  are  known 
except  In  the  first  case,  we  estimate  ^  =  *$(x).  The  true  (transformed)  coefficient  is  given  by 
D(y)  =  e  v.  In  Figure  1,  we  have  plotted  the  result  of  the  estimation,  using  M=5  for  the  parameter 
estimation  and  N  =6  for  the  state  approximation.  For  the  second  example,  we  estimate  =  ^(t),  with 
the  true  coefficient  given  by  =  D(t)  =  cos(t).  The  results  of  this  estimation,  again  using  M=5  and 
N=6,  can  be  seen  in  Figure  2. 

We  also  refer  the  interested  reader  to  [7],  where  an  example  of  the  Stefan  problem  appears.  As 
mentioned  above,  the  Stefan  problem  violates  (H2),  however,  the  numerical  method  is  nonetheless 
successful. 
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FIGURE  1.  Estimotion  of  D(y),  M=5 


Estimation  of  D(t),  M*5 
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